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I The steady state motion of cylindrical droplets under the action of external body force is 

, investigated both theoretically and via lattice Boltzmann simulation. As long as the shape- 

. invariance of droplet is maintained, the droplet's center-of-mass velocity linearly scales with 
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both the force density and the square of droplet radius. However, a non-linear behavior 
appears as the droplet deformation becomes significant. This deformation is associated 
with the drop elongation occurring at sufficiently high external forcing. Yet, independent of 
either the force density or the droplet size, the center-of-mass velocity is found to be linear 
in terms of the inverse of dynamic viscosity. In addition, it is shown that the energy is 



' mainly dissipated in a region near the substrate particularly close to the three phase contact 

- ' line. The total viscous dissipation is found to be proportional to both the square of force 



density and the inverse of dynamic viscosity. Moreover, the dependence of the center-of-mass 
velocity on the equilibrium contact angle is investigated. A simple analytic model is provided 
reproducing the observed behavior. 
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^ ; I. INTRODUCTION 

Individual droplets play a key role in many biological systems [l| . Droplet behavior is also crucial 
for numerous industrial applications such as in automobile manufacturing and drug production as 
well as glass industry. Consequently, understanding the underling physics behind droplet behavior 
and finding novel applications are currently an active field of research 0]. Recently, study 
of microdroplets has received lots of attentions both experimentally and by numerical modeling. 
For example, droplet spreading on chemically and topographically patterned substrates, droplet 
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evaporation, and wetting properties of superhydrophobic surfaces have been extensively studied 



in the hterature [4l4lll|. Particularly, controlling droplet motion is essential for many industrial 
purposes ranging from microfluidic devices to fuel cells and Inkjet printing jl^ . 



The equilibrium contact angle of a droplet, O^q, placed on a perfectly flat and homogeneous solid 
substrate is given by the Young equation, cos^eq = {(^sv — '^sl)/clv, where ar.v, crsr,, and asv 
are the surface tensions of liquid- vapor, solid-liquid and solid- vapor, respectively [l3]- However, 
in the case of moving drops, the advancing contact angle is often found to be larger than the 
receding IJ]. This deference can be considered as a measure of droplet deformation and it may 
appear in characterizing _of_ droplet velocity [isf] . A droplet may move due to a wettability or 



temperature gradient 



16l4l8l]. Recently, we were successful to report a spontaneous droplet motion 



191 ]. Obviously, a 



on a substrate topographically patterned with a step-wise gradient of pillars 
droplet may also move under the action of a body force, e.g., a falling drop on an inclined surface 
under the gravitational forcing. Depending on the material parameters of the considered system 
such as r], dynamic viscosity, and ^eq as well as superhydrophobicity of the substrate, droplets 



perform a sliding, rolling, or tank treading motion or a combination thereof [2CH22|| . Associated 
to a very high external forcing (i.e. sufficiently large velocity), drops may highly be elongated 
(pearling) and, further up, they exhibit a cuspid tail that emits smaller drops [23(. Introducing 
slippage at solid boundary is another issue that helps to characterize droplet motion 

Here, we concentrate on the steady state motion of cylindrical drops. However, despite the 
apparent simplicity of the problem, several issues, such as dependence of center-of-mass velocity, 
Ucm, and the dissipation loss on the material parameters and external forcing as fully as the role 



of droplet deformation are still not well understood |2j]. The steady state is reached due to the 



balance between the rate at which energy is imparted onto the droplet and the rate of energy 
dissipation. In general, there are different possible mechanisms for energy dissipation within a 
moving droplet: the viscous dissipation due to the velocity gradients, dissipation at the vicinity 
of the three phase contact line, and the dissipation in the precursor film which may form around 
the droplet in contact with a solid . Since, the numerical model used in the present studies 

does not take account of precursor film, we will focus on the effects related to dissipation only. 
This includes both the bulk of the drop as well as the vicinity of the contact line. Interestingly, as 
long as external force is sufficiently weak or -equivalently- the droplet volume is sufficiently low 
so that the drop approximately maintains its equilibrium shape during motion, the dependence 
of the center-of-mass velocity on external force and on the droplet volume can be easily worked 
out via simple rescaling of the relevant parameters. In particular, we find that the steady state 
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drop velocity is directly proportional to gR'^g/r], where g is the external body force (equivalent of 
the gravitational acceleration), iieff the effective drop radius and r] the shear viscosity. Deviation 
from this simple behavior is observed in the case of strong droplet deformation. However, since 
dynamic viscosity does not affect the droplet shape, the drop velocity remains proportional to 
1/r/ even in the strongly deformed limit. Using numerical simulations, we also calculate the local 
energy dissipation inside the droplet. It is observed that the main dissipation takes place within 
a volume below the drop's center-of-mass. Based on this observation, we propose a simple model 
which successfully captures the dependence of drop velocity on equilibrium contact angle. 



II. NUMERICAL MODEL 

Because of complicated nature of fluid flows, tractable analytical approaches are often limited 
to simplified systems. In addition, experimental studies are available only for a restricted range of 
parameters. In this context, computer simulations can help to bridge the gap between analytical 



approaches and experiments. In the past two decades, the lattice Boltzmann (LB) method 251432 1 
has proved itself as a powerful Navier-Stokes solver for simulating a wide range of complex fluidic 
systems. 

We employ a free-energy-based two-phase lattice Boltzmann (LB) model to solve the discrete 
Boltzmann equation (DBE) for the van der Walls fluid with the BGK approximation. A detailed 



description of the model can be found in references {33I . 



3^ . For the sake of completeness, however. 



a brief overview of the model is provided here. The DBE with external force F can be written as 

^I^A_f, .Vf _ fa- fa I (ea-u)-F 

dt v^"- A + pel •'^ ■ 

In the above, /a, Gq, and u are particle distribution function, the microscopic particle velocity 
and the macroscopic velocity, respectively. The parameter p stands for the fluid density, A is the 
relaxation time and Cs denotes the sound speed. The non dimensional relaxation time r = \/5t is 
related to kinematic viscosity hy v = Tc^6t. The equilibrium distribution function, fa'^, is given by 



f'a' = ^aP 



Gq, • U (Gq • u) u • u 



1 + ^ + 



(2) 



where WqAs a weighing factor. In order to eliminate the parasitic currents, the averaged external 
force experienced by each particle F is chosen in the potential form 



F = Vpc2-pV(/xo-KVV), 



(3) 
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where //q is the chemical potential and k the gradient parameter. The equilibrium properties of the 
present model can be obtained from a free-energy functional consisting of a volume and a surface 
part, 

^ = {e^{p) + ||V|2) dV - j{ct>ips)dS, (4) 

where V is the system volume and S the surface area of the substrate. The bulk energy density, Eq, 
can be approximated by Eq{p) = /3(/0 — /Ov)^((p — Pl)^) in which /? is a constant and both pL and py 
are saturation densities in liquid and vapor phase, respectively. The gradient parameter and the 
liquid-vapor surface tension can be computed as k = (3D'^{pi^ — py)'^/8 and a = {pi, — /3v)^\/2k/3/6, 
respectively. The interface thickness D, f3, pi^, and py are input parameters. The second integral in 
Eq. ([3]) is the contribution of solid-liquid interfaces in the total free energy ^. At equilibrium, there 
are two solutions that satisfy (pi = ib-y/2KE'o(/o). Minimizing the free energy functional ^ leads to 
an equilibrium boundary condition for the spatial derivative of fluid density in the direction normal 
to the substrate dJ-p = —(pi/K. The parameter (pi is related to 9eq via 

01 =^(PL-Pv)'sgn(f -0eq) 

x{cos(t)[l-cos(f)]}^/\ (5) 

where a = arccos(sin0eq)^. 




Figure 1: (color on-line) Difference between rescaled velocity fields, U2{y,z) — ui{y,z), for two different 
values of g (left), rj (middle) and i?cff (right) as indicated. Other control parameters of the simulation are 
as follows: — 0.16 and Rcs — 19.7 (left panel); g = 10^^ and i?cff — 19.7 (middle panel) and finally 
g = 5 X 10^^ and = 0.16 (right panel). The difference between rescaled velocity fields in computed after 
a shift operation such that the centcr-of-mass of the droplets coincide with one another. 

The advantage of this model for the current study is both the possibility of achieving a high 
density ratio and, as it was already mentioned, the elimination of parasitic currents at the liquid- 
vapor interface. It is important to note that the elimination of the spurious currents is an important 
step towards a reliable description of fluid dynamics inside a droplet. Simulating a two-phase system 
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with a high density ratio, on the other hand, not only is more reahstic but also allows to significantly 
reduce the finite size effect related to the dissipation loss in the vapor phase. 

In our LB simulations, the bounce-back rule is imposed at solid boundaries. For the open 
boundaries (in the x and z-directions), the periodic boundary condition is applied. A body force, 
pg, is applied to the liquid phase along the z-direction. The body force, however, monotonously 
decreases through the interface and it vanishes in the gas phase. This accoTints for the fact that 
the gas remains inert (static equilibrium) in the limit of zero droplet size. All the quantities in this 
paper are given in dimensionless LB units. The parameter the interface thickness D and the 
saturation densities are fixed to 0.01, 5, 1, and 0.01, respectively. This choice of the parameters 
leads to a surface free energy of cr 2± 0.004. Note that, in order to focus on situations, which can 
be easily controlled in real experiments, we do not change surface tension or liquid density in our 
simulations. Depending to the case of interest, the parameters r, ^^cq, RcS-, ^-nd 9 he in the ranges 
[0.02,1.6], [35°, 140°], [22,75] and [10~^, 10"^] in the order given. Typically, we use a simulation 
box of size Lx x Ly x = 2 x 120 x 120 lattice nodes. However, for large droplets, we increase the 
size of the simulation box (in the y and ^-directions) ensuring that there are no finite size effects. 
The volume of droplet is given hj V = SL^ where S is the surface of droplet's cross-section normal 
to the x-direction. For the cylindrical geometry considered in this study, we define the droplet's 
effective radius as Res = (S/tt)^/^. 

III. A SIMPLE SCALING RELATION 

Here, we investigate the effect of external forcing on the steady state motion of cylindrical drops 
on a fiat surface. In addition, the infiuence of system parameters such as droplet size ReS, viscosity 
r) as well as equilibrium contact angle ^eq on the steady state velocity of the droplet's center-of-mass 
is addressed. 

By driving a droplet via an external body force, we mimic a real situation in which a droplet 
moves downward on an inclined surface due to the gravity. The external force does work on droplet 
with a rate equal to the total force applied on the droplet multiplied by the droplet's center-of-mass 
velocity, gpVUcm- In the steady state, this energy is entirely transferred into dissipation. On the 
other hand, the total viscous dissipation is given by Jy SijaijdV = Jy aijaij / {2i])dV , where the 
strain rate and stress tensors are given by Sij = {dui/dxj + duj/dxi)/2 and aij = 2r)Sij, the later 
relation being valid for non-diagonal (shear) components of a Newtonian fiuid (note that, due to 
the incompressibility of the liquid phase, the diagonal components of the strain rate and stress 
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tensors are not relevant here). One thus obtains 



gpVUcm = ^ y (^ij(^ijdV = '^'n j SijSijdV. (6) 

As long as the shape of the droplet does not change, it is reasonable to take R^s as a characteristic 
length. We also chose C/cm as a characteristic velocity and introduce dimensionless quantities such 
as X = Xa/Rcs and Ua = Ua/Ucm- Using these rescaled quantities, the strain rate tensor can also 
be written as Sij = Ucia/ Rcsidui/dxj + duj / dxi) /2 = Ucm/ ReS^ij- Inserting this into Eq. ([6]) 
yields 

gpVUcm = -^S^ / S,,S,,dV. (7) 

where we also made the volume element dimensionless {dV = VdV). The important step is now to 
assume that the rescaled velocity field within the droplet does not change upon a variation of the 
external force, drop radius or viscosity provided that the shape of the droplet remains constant. 
With this assumption, the integral in Eq. ([7]) becomes a constant 'shape factor' and one obtains 

f/c. oc (8) 
V 

It is noteworthy that, in the above model, the dependence of Ucm on R^g arises from the rescaling 
of the strain rate tensor Sfj only. In particular, it remains valid regardless of the dimensionality of 
the space. Interestingly, when expressed in terms of droplet volume, V, the spatial dimension, d, 
does play a role. This is simply a consequence of the fact that Rcs oc V^^'^. In particular, Ucm oc V 
in 2D, while Ucm oc V^/^ in 3D. 

In order to test the above assumption of the scale invariance, we have performed a series of 
lattice Boltzmann simulations while varying g, r] and Rcs in a range where droplet shape remains 
unchanged. The simulated velocity fields are then compared with one another by first rescaling 
the relevant velocity and length scales (see the text below Eq. ([6])) and then plotting the difference 
of the thus obtained velocity fields. 

Results of such an analysis are illustrated in Fig. [TJ As seen from this figure, the rescaled velocity 
fields are very close to each other almost in the entire droplet with deviations in the vicinity of 
the three phase contact line. Noting that these deviations (being at most of the order of 10%) 
are limited to a small fraction of the droplet's volume, the relative contribution of these deviations 
to the integral in the right hand side of Eq. ([7]) becomes quite negligible in all the cases shown. 
Obviously, the assumption of a scale invariant velocity field is a good approximation to the actual 
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Figure 2: (Color on-line) [/cm versus body force g for different values of the effective droplet radius RcS as 
specified. A linear behavior is visible at sufficiently low g. The range of the validity of this linear regime is 
progressively restricted as droplet radius increases. In the right panel, jyC/cm/ (ffP^off) plotted versus g for 
exactly the same data as in the left panel. In all these simulations, shear viscosity and equilibrium contact 
angle are set to 77 = 0.16 and 9y = 90°. 
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Figure 3: (Color on-line) Ucm versus Rl^g for different choices of the body force g as indicated. Again, a 
linear behavior is visible at sufficiently low i?cff- The range of the validity of the linear behavior shrinks 
upon a raise of the body force. Following the same idea as in the right panel of Fig. [51 we plot in the right 
panel ijUcm/ (gpRcs) versus R^g for exactly the same data as in the left panel. In all these simulations, shear 
viscosity and equilibrium contact angle are set to r/ = 0.16 and Oy = 90°. 



flow behavior in the studied parameter range. Equation ([8]) is thus expected to well describe our 
data as long as droplet shape is unaltered. 

The presence of a parameter range for the validity of Eq. ([8|) is evidenced in Fig. [2l where the 
center-of-mass velocity, Ucm, is depicted versus force density, g, for different droplet sizes. As seen 
from this figure, the range of the validity of scaling relation Eq. ([8|) extends to larger g as droplet 
size decreases. Conversely, the larger the droplet, the earlier the onset of significant deviations. A 
similar trend is also observed in Fig. [3l where droplet size is varied as control parameter for three 
different choices of g. 

In the present study, the shape of droplet is determined by the competition between the sur- 
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face force and the total body force. For a cylindrical drop of the cross sectional radius i?efr and 
axial length L^, this leads to au/L^ < gpR^^L^ as a condition for a significant deformation. By 
introducing the Bond number, Bo = pgR'^g/o'i^^, one sees that strong deformation is expected for 
Bo > 1. Within prefactors of the order of unity, the same condition for drop deformation is also 
obtained in the case of a spherical droplet (to see this, replace by -Rcff)- It must be emphasized 
here, that this criterion is based on a scaling argument and the precise value of the Bond number for 
the transition from undeformed to a deformed state may be different from unity. What is essential 
here is the fact that a higher Bond number leads to a higher degree of deformation. In the case of 
our simulations, for example, slight but observable deformation occurs already for a Bond number 
as low as 0.25 (Fig. [Hb) with a significant increase in the deformation state as Bo increases from 
0.25 to 0.72 (Fig.Ht). 

Droplet shapes and the corresponding momentum fields are shown in Fig. H] for three typical 
values of g. As seen from the left panel of this figure, for a sufficiently weak body force (here 
g = 10~^), the deformation of the droplet is quite negligible but it becomes important upon an 
increase of g (middle and right panels). 
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Figure 4: Droplet shape and the corresponding momentum field in the center-of-mass frame for three different 
values of the driving force g. As g increases, the deformation becomes more pronounced. In the left panel, the 
deformation is negligible and the droplet 's center-of-mass velocity obeys the simple relation Eq. ([8]) for driving 
forces below the specified value. The middle panel marks the onset of deviations from Eq. (jH) and the left 
panel is well beyond the validity of this simple scaling relation. A rolling motion is clearly visible regardless 
of the deformation state of droplet. In all the cases shown, the droplet's effective radius, dynamic viscosity 
and equilibrium contact angle are fixed to i?cff — 19.7, rj = 0.16 and 0y = 90°, respectively. Recalling that 
c^LV = 0.004 and pL = 1-0, the Bond number from left to right reads Bo ~ gpR^^/ aij\r « 1 x 10~^, 0.25 
and 0.72. 



Furthermore, Fig. d] also shows the momentum field inside the droplet providing direct evidence 
for the existence of rolling motion in the center-of-mass frame of reference. Thus, an observer 
moving with the droplet's center-of-mass will confirm the presence of a well established rolling 
motion inside the droplet regardless of its deformation state. This rolling motion is associated to 
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Figure 5: (Color on-line) Ucm versus 1/t] for droplets with different degrees of deformation. The labeles (a)- 
(c) refer to deformation states shown in Fig. |4l respectively (the velocity in the case (a) has been multiplied 
by a factor of 10 for better visibility). In all the cases shown, a perfect linear variation is seen in accordance 
with Eq. ([5]) . In the left panel, we plot simulation results for two different choice of i?off and g but keeping 
the product Rlf^g almost unchanged. In this case, the velocities of both droplets fall onto a single line which 
also confirms the validity of Eq. ([5]). d^q is fixed to 90° in all cases. 
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22l |. Interestingly, similar rolling 



the tendency of droplet to minimize its total dissipation loss 

motion are also observed in molecular dynamics simulations of polymeric liquids |24 

We close this section by addressing the effect of viscosity on Ucm- For this purpose, we mention 
that a change in viscosity only affects the time scale of the entire simulation. In particular, a 
variation of viscosity has no influence on the shape of droplet. Consequently, we expect Ucm oc l/rj 
regardless of the deformation state of droplet. This expectation is confirmed in Fig. [5l where Ucm 
versus 1/r/ is shown for droplets with different degrees of deformation. 



IV. LOCAL VISCOUS DISSIPATION 

In this section, we provide a detailed analysis of the local dissipation rate, </>(t') = a'^{r)/2r] 
inside droplet. As will be shown hereafter, the insight gained via these investigations enables us 
to propose a simple model capable of accounting for the dependence of the total dissipation rate, 
$T = / 4>{r)d^''r, on the contact angle, 6y- Equating this to the work done by the external force 
then yields a relation between the droplet 's center-of-mass velocity and the equilibrium contact 
angle. 

It is noteworthy that, unlike conventional Navier-Stokes solvers, the lattice Boltzmann method 
does not require — although allows for — the computation of velocity gradients to obtain the local 
stress tensor. Rather, it offers the unique possibility of obtaining the stress tensor locally via the 
non-equilibrium part of the populations. In this regard, particular attention has been payed to a 
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correct implementation of the stress computation [38 1 . 

In order to figure out at which parts of droplet the energy is mainly dissipated, we compute 
local dissipation rate along the three lines labeled by A, B, and C in the panel (a) of Fig. [6l The 
variation of (p along these lines is depicted in the next panels of Fig. [6l We first note that (p 
is negligible in the gas phase, which is often the case in real experiments due to the low vapor 
pressure. Furthermore — as a comparison of the panels (a), (b) and (c) reveals — the strongest 
dissipation occurs in the vicinity of the three phase contact line (panel (c)), which is roughly two 
orders of magnitude larger than the dissipation rate inside droplet (panel (b)). This behavior can 

be rationalized due to the fact that large velocity gradients occur near the substrate particularly 

1 i 

in the vicinity of the three phase contact line ||39]. However, one must realize that bulk dissipation 
acts in a larger domain than the dissipation close to the contact line and thus may eventually 
dominate the overall dissipation rate if the droplet is sufficiently large. 

An interesting feature, relevant for our further analysis is the fact that viscous dissipation 
inside droplet is mainly localized to regions below the droplet 's center-of-mass (panel (a) in Fig. [6]). 
This idea is further evidenced in Fig. [7] (a), where we plot the viscous dissipation integrated 
along a horizontal line, ^{y) = j^'^ 4>{y,z)dz as a function of vertical position y (distance from 
the substrate). Indeed, as expected, viscous dissipation mainly occurs in a region specified by 
y < Ycm- This finding is further underlined by showing in Fig. [7] (b) the relative contribution, 
^r(2/)) to total dissipation within a region restricted between the substrate and a horizontal line 
at y, ^R{y) = f^^{y')dy'/ j^^ <^{y')dy'). It is visible from Fig.[7](b) that 96% of total dissipation 
occurs in a region specified by y < Y^m- 



H 


A 


B 


^ 1 








1 C 



(a) 



20 40 60 80 100 120 
z 



I 

_~ 0.8 
2 0.6 

X 

-0.4 
0.2 




along the 
line A 




(b) 



15 30 45 60 

y 




5 

, 4 

3 
' 2 

1 

oL . . . „ 

30 60 90 120 



along the 
lineC 



(d) 



Figure 6: (a) Illustration of the typical shape of a droplet and the lines along which viscous dissipation 
is determined ('CM' stands for the center-of-mass). The system parameters are g = 10^^, i?cff = 28.2, 
T] = 0.16, and = 90°. (b-d) The variation of local viscous dissipation rate, (j) — (j"^ [r) / {2rf) , along lines 
A, B and C as indicated. Note that is almost negligible in the gas phase. It has a large value close to the 
substrate, but rapidly decreases far from the substrate. It is also highly enlarged in the vicinity of triple 
contact line. 



Before working out an important consequence of this observation, we first check whether it 
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Figure 7: (Color on-line) variation of the local dissipation rate $ and the relative dissipation rate as 
a function of y, corresponding to the droplet shown in Fig. [6] in the left and right, respectively. The main 
contribution to the total dissipation occurs close to the substrate in a region given by y < ycm- 



remains valid upon a variation of shear viscosity and driving force. Inserting Eq. 
hand side of Eq. ([7]) , it is seen that the total dissipation rate is expected to obey 



9„2 „2 p2+d 



in the right 



(9) 



In order to verify Eq. Q, we determine $(y) for different values of body force and viscosity. 
Typical plots of the thus obtained results are shown in Figs. [8] and [H These data clearly underline 
the validity of Eq. ([9]) within the studied range of parameters. 
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Figure 8: Left: Dissipation rate integrated along a horizontal line at y, Right: The same quantity as 

in the left panel but divided by g^. The observed master curve supports the validity of Eq. 



In addition to supporting the validity of Eq. Q, the data shown in Figs. [8] and [9] provide 
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Figure 9: A similar plot as in Fig. [5] but now for various fluid viscosities rj. Here, the right panel depicts 
<i>(2/)-data from the left panel multiplied by rj. Again, the validity of Eq. ([9]) is supported by the master 
curve. 



further evidence for the fact that most part of dissipation occurs in the region below the droplet 's 
center-of-mass. Based on this observation, we propose a simple relation allowing to describe the 
dependence of droplet velocity on contact angle. 

Our simple analytic model is based on scaling arguments. To proceed, we start with the 
energy balance equation for a cylindrical droplet of axial length in the steady state. Us- 
ing the translation invariance with respect to the x-coordinate, one can write gpirR'^gLxUcm = 
Lx{ri/2) Jq j 'j'^dzdy, where 7 is the local shear rate. Since the energy is almost completely dissi- 
pated in a region below Ycm, we can safely restrict the upper limit of the integration to Ycm and 
rewrite the energy ballance equation as 2gpTTR'^gUcm = V /q^"" / "j'^dzdy (see Fig. [T]). Neglecting 
droplet deformation, the droplet 's cross-section is a circular segment with a base contact angle of 
6eq- Here, we assume that 7 simply scales as C/cm/^cm throughout the droplet. This may appear 
as a crude approximation, but it allows to obtain a solvable analytic expression. Furthermore, we 
approximate the surface of the droplet below as that of a rectangle of height and length 
Iz- Adopting this, the right hand side of the energy balance equation can now be estimated by 
{r]/2)lzYcm X (Ucm/Ycmf = {r} / 2)1^11"^ J Y^m- Thus, one obtains 2gpTTRlgUcm = v^zU^m/^cm, which 
then yields Ucm = gp'^Rcs^cm/i^zV)- the considered geometry, the quantity Ycm/lz, is only a 
function of ^oq- Taking this into account, we finally arrive at 



c 



9PR 



eff 



4sin^6', 



cq 



3(2eeq - Sin2eeq) 



cotO. 



eq 



(10) 



The validity of the model has been tested in Fig. [10] for two different droplet radii. This simple 
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model reproduces well the simulation results. Interestingly, the fitting prefactor, C, for both 
investigated droplet sizes are very close to each other (0.56 and 0.57) showing the consistency of 
the model. 

We would like to emphasize that the present approach is different from conventional approaches, 
where the integration is taken over the entire volume of droplet. Following the conventional route, 
one would obtain a different expression, Ucm = C'(5'/5^cfT/^)(-'^~*^o^^cq)^/(^cq — sin^cqcos^eq), which, 
as shown in Fig. [TOl is not successful in capturing the observed behavior. 

It is noteworthy that an extension of Eq. (jlOp to 3D can simply be obtained by writing the 
energy ballance equation in 3D and replacing the corresponding expression for Ycm/lz by that of a 
spherical cap. 
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Figure 10: Droplet velocity versus equilibrium contact angle for two different droplet volumes as indicated. 
Full solid lines are best fit results to Eq. ([TOl) while dashed lines give best fit results to Ucm — C{gpRlg/ri){l — 
cos0oq)^/(^^cq — sin(?cqCOS0cq) (sce the text). The force density and dynamic viscosity are fixed to g = 10^"^ 
and 77 = 0.16 respectively for both droplets. 



V. CONCLUSION 

We use a two-phase lattice Boltzmann method to study the dynamics of cylindrical droplets 
on a flat substrate under the action a gravity-like external force density. Starting from the energy 
ballance equation, we first drive a simple analytic relation, Eq. ([8]), indicating that — as long as 
the shape- invariance of droplet is maintained — droplet 's center-of-mass velocity, linearly scales 
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with force density, and the square of the droplet radius. At strong body forces or large droplet 
volumes, deviations from Eq. ([8]) are observed. A survey of droplet shape within our simulations 
suggest that droplet deformation is indeed the main cause of observed deviations from the simple 
scaling relation. Interestingly, however, the droplet's center-of-mass velocity remains proportional 
to the inverse of the dynamic viscosity regardless of droplet's deformation state. This is in line 
with the idea that viscosity merely affects the time scale of the problem with no influence on 
droplet shape. A detailed study of the local dissipation inside droplet is also provided. A results 
of these investigations is that dissipation mainly occurs close to the three phase contact line and 
within a region below the droplet's center-of-mass. Using the latter observation, we propose a 
simple analytic expression accounting for the dependence of droplet velocity on the equilibrium 
contact angle. Results of computer simulations confirm the validity of this simple model. 
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